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S ■ Abstract 

O ■ We study theoretically Electron Paramagentic Resonance (EPR) linewidths 

for single crystals of the molecular magnets Feg and Mni2 as functions of 
energy eigenstates Mg, frequency, and temperature when a magnetic field 
along the easy axis is swept at fixed excitation frequency. This work was 
motivated by recent EPR experiments. To calculate the linewidths, we use 
density-matrix equations, including dipolar interactions and distributions of 
the uniaxial anisotropy parameter D and the Lande g factor. Our calculated 
linewidths agree well with the experimental data. We also examine the line- 
shapes of the EPR spectra due to local rotations of the magnetic anisotropy 
axes caused by defects in samples. Our preliminary results predict that this 
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effect leads to asymmetry in the EPR spectra. 
PACS numbers:75.45.+j,75.50.Xx,76.30.-v 
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Molecular magnets such as Mni2^ and Fes^ have recently drawn much attention because of 
the macroscopic quantum tunneling of their magnetizations at low temperatures^ and their 
possible applications to quantum computing.'' These materials consist of many identical 
clusters (see Fig. 1) with the same magnetic properties and characteristic energies. Each 
cluster has many different species of ions and atoms, with a total spin angular momentum 
in the ground state of 5=10. The clusters have strong crystal-field anisotropy and, thus, a 
well-defined easy axis, and the magnetic interaction between different clusters is weak. 

Recently, multi-frequency Electron Paramagnetic Resonance (EPR) measurements^"^ 
on single-crystals of the molecular magnets Feg and Mni2 showed interesting results in the 
linewidths as functions of the energy level Ms, excitation frequency, and temperature when a 
magnetic field along the easy axis was swept with the excitation frequency fixed. Figures 2(a) 
and (b) show the experimental linewidths as a function of the value of Mg from which the 
spin system is excited. For example, Ms—5 in Fig. 2(a) denotes the transition Ms=5— >4. 
The experimental results are that (i) the hnewidths are about 200 G to 1400 G at T=10 
K for Feg [Figs. 2(a) and (b)], and about 1000 G to 2000 G at T=25 K for Mn^] (ii) the 
linewidths increase non-linearly as a function of the absolute value of the energy eigenstate 
Ms [Figs. 2(a) and (b)]; (iii) the hnewidths attain a minimum at Mg—l and Mg—O [Fig. 2(a)]; 
(iv) the linewidths decrease with increasing frequency (compare the linewidth for Mg=5 in 
Fig. 2(a) with that in Fig. 2(b)). 

In this paper, we present the theoretical EPR hnewidths for single crystals of Fcg and 
Mni2, using density- matrix equations with the assumption that the uniaxial crystal- field 
anisotropy parameter D and the Lande g factor are randomly distributed around their mean 
values ("£)-strain" and "g'-strain" effects^) due to possible random defects and impurities in 
the samples. We find that the calculated linewidths^ agree well with recent experimental 
data^'^ as functions of the energy level Ms and the frequency. We also briefly consider the 
effect of the temperature on the linewidths. Finally we present a preliminary analysis of 
how local rotations of the magnetic anisotropy axes, caused by defects in the samples, affect 
the lineshapes. 
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To obtain the linewidths for Feg, we consider a single-spin system with 5"= 10 in a weak 
oscillating transverse field. We choose the easy axis to be the z-axis. Since the Feg clusters 
have an approximate D2 symmetry, the lowest-order ground-state single-spin Hamiltonian 



Ho = -DSl - E{Sl - Si) - qixbH^S, , (1) 

where D=0.288kB and the transverse anisotropy parameter, E=0.0A3kB-^ Here Sa is the a- 
th component of the spin operator, g is the Lande g^-factor (~2), //^ is the Bohr magneton, 
and Hz is the longitudinal static apphed field. When is large enough, the eigenvalue of S^, 
Mg, is a good quantum number. Next we introduce an interaction V{t)=VoCOs{ut) between 
the spin system and an oscillating transverse field with angular frequency a;=27r/, where 
VqocHxSx. We treat V{t) as a small perturbation to Hq. The interaction between the spin 
system and the surrounding environment can be understood by density-matrix equations. 

We consider the case that the frequency cu is fixed while is quasi-statically varied to 
induce a resonance. With the selection rule AM^^il, solving the density-matrix equations 
for the population change with time in the level Mg due to V{t), to first order in Vq near 
resonance, provides the power absorption between the levels Mj and Mg — 1:^^ 



^ V A(pm„m,-Pm,-i,m,-i), (2) 
h 

A=. 



dt 

-1,M. 



Here Em, is the energy of the level Mg, V=\{Mg\Vo\Ms — p{t) is the density matrix of 
the spin system. The subscripts represent eigenstates of the longitudinal part of Ho, p{t)m'm 
= {m'\p{t)\m), 'ym',m={Wm + Wm')/2, Wm=J2k^m Wkm, and Wkm is the transition rate^^ from 
the m-th to the k-th eigenstate, i^res is the resonant field, and Ti^Ms-i,Ms/ QPb gives a natural 
hnewidth which is about 5 to 50 G at 10 K and increases as Mg decreases. However, the 
experimentally observed linewidths are much larger than the natural linewidths and decrease 
as Mg decreases [Figs. 2(a) and (b)]. 



To resolve these large discrepancies, we first assume that D and g are independent random 
variables with Gaussian distributions centered at 0.288/cb and 2.00, with standard deviations 
ud and cr^, respectively. Then we calculate the average power absorption at a fixed frequency 
and 7=10 K by integrating over Eq. (2) using Mathematical^ to obtain the hnewidth as a 
function of Mg. The lineshapes depend on the magnitudes oi (Td and Ug, compared to the 
natural linewidth determined by 7Ms+i,Ms- If ctd and Ug are much larger than (comparable to) 
the natural linewidth, then the lineshape is Gaussian (Lorentzian). At intermediate values 
of (Td and Ug, the absorption lineshapes are neither Gaussian nor Lorentzian. Our numerical 
calculations show that the distribution in D narrows the linewidths linearly with decreasing 
absolute value of 2Ms— 1 [Figs. 2(c) and (d)], with a slight rounding close to the hnewidth 
minimum (Ms=l and M^^O) [Fig. 2(c)]. On the other hand, the distribution in g broadens 
the linewidths with decreasing Mg because the resonant field increases with decreasing Ms 
[see Eq. (3)].^'^ For small- M, transitions, the hneshapes are close to a Lorentzian because 
the natural linewidths are not very small compared to the measured linewidths. 

Next, we consider the effect of the dipolar interactions between different clusters. (There 
is no distribution in cluster size.) The dipolar interactions narrow the linewidths as M, 
decreases because the resonant field becomes stronger for smaller- Mg transitions, and the 
stronger resonant fields lead to a more polarized system [Figs. 2(c) and (d)]. The details of 
this effect on the linewidths at a fixed temperature were reported elsewhere.^ The dipolar 
interactions give rise to a temperature dependence of the linewidths and a shift in the 
positions of the resonance lines, because at low temperatures {kBT<^hui) high energy levels 
are not populated. Our study shows that a linewidth for a particular transition increases 
and then smoothly decreases with increasing temperature. The maximum linewidth as a 
function of temperature moves towards lower temperatures for larger- transitions, which 
agrees with experiments. 

The competition between £)-strain, gr-strain, and dipolar interactions determines the 
overall features of the linewidth, as a function of Mg. We have varied ctd, ag, and the 
effective distance between neighboring dipoles, d, within an experimentally acceptable range 
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in order to fit the experimental data. 

For the Feg sample examined, the calculated linewidths agree well with the experimental 
data at the measured frequencies (/ = 68, 89, 109, 113, 133, and 141 GHz) at T=10 K, using 
(7d~0.01£) and d^l2 A. [Figs. 2(a) and (b)] As shown in Figs. 2(c) and (d), the D-stredn 
effect and the dipolar interactions are equally important for the linewidths of the sample, 
while the gi-strain does not contribute significantly (not shown). 

For Mni2, we perform a similar analysis as above, using the ground-state single-spin 
Hamiltonian^° 



with D^0.55kB, C=1.17 x 10~^kB, and ^=1.94.^^ Detailed analysis can be found in the 
literature.^ For this sample, the D-strain and gi-strain effects play significant roles in the 
linewidths, while the dipolar interactions are not as important as for the Fes sample. 

Next we discuss how the distribution of the directions of the magnetic anisotropy axes 
of clusters affect the EPR lineshapes. Because of defects in the samples, each cluster sees 
a slightly different crystal field due to the surrounding clusters, compared to the situation 
in a perfect crystal. We assume that these shghtly different crystal fields result in local 
rotations of the magnetic anisotropy axes of some of the clusters by an angle 6 from the 
crystal c axis. The angle 9 is assumed to have a Gaussian distribution about zero with a 
small standard deviation. Hereafter a, b, and c denote the crystal axes, while x, y, and z 
denote the magnetic anisotropy axes of a single cluster. 

As a preliminary study, we examine the lineshapes for Mni2 when the magnetic field is 
applied along the c axis and the magnetic anisotropy easy axis of a single cluster is tilted by 
6 from the c axis. Then the single-spin Hamiltonian, in terms of the spin operators along 
the magnetic anisotrpy axes, becomes 



For simplicity, we drop the fourth-order anisotropy terms. We also set -0 = 90°, where i/j is 



Ho^-DSl-CSl-giXBH^S, 



(4) 



H — —DS^ — q/ibHSz cos 9 + q/ibHSx sin 9 . 



(5) 
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the third of the Euler angles defined in hterature,^^ which does not affect the eigenvalues of 

n. 

Because of the local rotations of the magnetic anisotropy axes, some of the clusters 
experience slightly different resonant fields than those aligned with the crystal c axis. In 
contrast to the Gaussian distributions of the resonant field due to the distributions in D and 
gi, this effect gives rise to an asymmetry in the resonant field, which leads to an asymmetry 
in the EPR hneshapes. Assuming that 9 is small, we treat the terms proportional to sin 9 as 
small perturbations to the rest of the terms in Eq. (5). Using second-order perturbation 
theory, we obtain the resonant fields as a function of 6*, shown in Fig. 3. Each curve 
in Fig. 3 is symmetric about ^=0. Figure 3 shows that for transitions between and 
Mg+l (Ms=— 4^— 3) an asymmetry appears in the direction of decreasing field, while for 
transitions between Mg and Mg—1 {Mg=3^2, Mg=2^1) the asymmetry is in the direction 
of increasing field. At the examined frequency (/=66.135 GHz) the asymmetry effect is more 
significant for the small- transitions. Based on the analytic form of the resonant field, we 
can obtain analytically the distribution of the resonant field at a particular frequency for 
transitions between the levels Mg and M^il, including Gaussian distributions of D, g, and 
9}^ For simplicity, we have neglected the effects of natural linewidths, dipolar interactions, 
and temperature on the EPR lineshapes. 

In conclusion, we have examined the EPR linewidths as functions of Mg, frequency, and 
temperature for single crystals of the molecular magnets Fcg and Mni2. We found that the 
distribution in D is important to explain the linewidths for both of the molecular magnets, 
and that the dipolar interactions are crucial to understand the temperature dependence of 
the linewidths. Our preliminary study of the EPR lineshapes shows that the distribution 
of the directions of the magnetic anisotropy axes of clusters provides an asymmetry in the 
EPR spectra. 

Funded by NSF Grants DMR-9871455, DMR-0120310, and DMR-0103290, and FSU- 

CSIT, FSU-MARTECH, and Research Corporation. 
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FIGURES 




FIG. 1. Schematic diagram of the magnetic core of (a) the Fes cluster and (b) the Mni2 cluster. 
Each Feg cluster has D2 symmetry, while each Mni2 cluster has tetragonal symmetry. Both clusters 
have ground-state spin S = 10 (for Fes, 5 = (6 - 2) x 5/2, and for Mni2, 5 = 8x2-4x 3/2). 
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FIG. 2. (a) and (b) Experimental (filled circles) and theoretical (open circles) Pull Width at 
Half Maximum (FWHM) vs Mg for different frequencies at T=10 K for Fes- (c) and (d) Line 
broadening due to the Z)-strain (filled squares) and dipolar interactions (open squares) for (a) and 
(b), respectively. Here H\\z. The fit is best when the standard deviation of D is O.OID and the 
effective distance between dipoles is 12 A. 
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FIG. 3. The resonant field vs 6 for the transitions (a) Ms=-A^-3 and Ms=3^2, and (b) 
Ms=2— >1 at /=66.135 GHz for || c. Each curve is symmetric about ^=0. 
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